THE STRUCTURE OF COLD DARK MATTER HALOS: RECENT 
INSIGHTS FROM HIGH RESOLUTION SIMULATIONS* 



MARCEL ZEMP 

University of Michigan, Department of Astronomy 
500 Church Street, Ann Arbor, MI 48109-1042, USA 

mzemp @umich. edu 

Received 2. September 2009 

We review results from recent high resolution cosmological structure formation simula- 
tions, namely the Via Lactea I & II and GHALO projects. These simulations study the 
formation of Milky Way sized objects within a cosmological framework. We discuss the 
general properties of cold dark matter halos at redshift 2 = and focus on new insights 
into the structure of halos we got due to the unprecedented high resolution in these 
simulations. 
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1. Introduction to Cosmological Structure Formation Simulations 

In the standard cosmological framework, the two most dominant components in the 
universe are dark energy (DE) and cold dark matter (CDM). The nature of dark 
energy is unknown and in the standard model it is assumed that it is a vacuum 
contribution in the form of a cosmological constant (generally denoted by A) in 
Einstein's field equation of general relativity Also the nature of the CDM is not 
understood and it could consist of weakly interacting massive particles (WIMPs) 
that only interact via the weak nuclear force and gravity. A possible WIMP candi- 
date is the lightest supersymmetric particle (LSP) that is predicted in supersym- 
metric models of particle physics .1212141 q \^ means th a t the thermal velocities of 
the CDM particles are small which leads to the formation of tiny structures in the 
universe.'^' The baryonic matter, of which all the stars and planets are made of, 
contributes only ca. 17% of the total matter content.^ This standard cosmological 
framework is also known as the concordance or ACDM model. 

The parameters of the ACDM model are known to a good precision.'^ This leads 
to well defined initial conditions for cosmological structure formation simulations, 
which are completely determined by these parameters. The initial conditions are 

"This invited review is based on talks given by the author in 2008-2009. 
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set up at early cosmological times when the density perturbations, which probably 
originate from quantum mechanical fluctuations in the inflationary epoch,^ are still 
in the linear regime. 

In most cosmological structure formation simulations, the baryonic matter is re- 
placed by CDM due to computational reasons. The baryonic physics is complex and 
difficult to model in a self consistent way, which makes more realistic simulations 
very expensive. But the replacement is physically well motivated since for the large 
scale clustering properties of the CDM halos the baryonic contribution can be ne- 
glected. The baryons are only important at the centers of CDM halos: baryons can 
dissipate energy through radiative processes, hence they concentrate at the centers 
of the potential wells formed by the CDM halos where they can lead to a significant 
alteration of the CDM distribution. 

The non-linear evolution of the CDM is modelled as a collisionless fluid that 
interacts via gravity in an expanding universe. Collisionless means that the CDM 
can be described by a continuous distribution function that obeys the collisionless 
Boltzmann equation. The collisionless Boltzmann equation is a consequence of mass 
conservation under a Hamiltonian flow. The gravitational potential is then calcu- 
lated via a Monte Carlo sampling of the continuous CDM distribution function by 
N particles/bodies (sampling points). Hence this technique is called the A^-body 
method.^ The sampling of the continuous distribution function introduces artificial 
singularities that are cut-off by softening the singular point particle potential at a 
softening scale e.^ 

In the last few years, an enormous effort in our collaboration went into under- 
standing the structure of CDM halos of the size of our own galaxy, the Milky Way. 
In several large simulation projects, we put all the effort in resolving one such ob- 
ject within a larger cosmological volume, with a ACDM cosmology, with as much 
particles as computationally possible. This is achieved via a multi-mass nested re- 
finement technique EEH The first simulation in this series was the Via Lactea I simu- 
l ation mMI5MmMMMl wh ich was followed by Via Lactea nMMnMMMMMMB 
Finally, the GHALO project was the first simulation to resolve a Milky Way sized 
CDM halo with more than one billion particles.^ The Aquarius project is a sim- 
ilar simulation project that was done by a different group.^* 32 * 33 * 34 ! Aquarius has 
an equivalent resolution to GHALO. In this brief review, we will mainly focus on 
results from our projects. 

2. Numerical Aspects 
2.1. Initial conditions 

The initial conditions of Via Lactea I & II and GHALO are based upon the Wilkin- 
son Microwave Anisotropy Probe (WMAP) three-year data^ and the three simula- 
tions have very similar cosmological parameters. In general, the initial conditions 
for cosmological structure formation simulations consist of a cubic volume of side 
length Lbox (see Table [T]) with periodic boundary conditions. 
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Creating initial conditions with billions of particles is a real challenge. In the 
case of GHALO, the GRAFIC1 and GRAFIC2 packages,^ which are often used by 
the community to produce initial conditions]^] had to be parallelized The initial 
conditions are then set up at high enough redshifts zic (see Table [T]) when the 
density perturbations are still in the linear regime. The softening lengths e, where 
the singular point particle potential is cut-off, is also given in Table [I] 

2.2. Time evolution 

The consequent time evolution is done by numerically solving the coupled collision- 
less Boltzmann and Poisson equations. In the case of Via Lactea I & II and GHALO 
this was done with the code PKDGRAV.E5ES1 The latest version of PKDGRAV is a 
tree code that uses a fast multipolc method (FMM) in Cartesian coordinates in or- 
der to calculate the gravity.lMSSI PKDGRAV uses a 5th-order reduced expansion in 
potential for faster and more accurate force calculation in parallel, and a multipolc 
based Ewald summation technique for periodic boundary conditions.^ Much bet- 
ter parallel computing efficiency was achieved by using Single Instruction Multiple 
Data (SIMD) vector processing. 

For the time integration, PKDGRAV uses a kick-drift-kick leapfrog integration 
scheme which is a second order, explicit, time-reversible and symplectic integrator 
when used with constant time-steps The individual time-steps of the particles 
are chosen according to their local dynamical time which leads to a more accurate 
orbit integration in the centers of CDM halos and an overall more optimal time-step 
distribution which as a consequence requires less force evaluations.^ 

3. CDM Halos 

3.1. General description 

In a ACDM cosmology, the main host halos contain a wealth of subhalos - smaller 
bound structures similar to their host halos. These subhalos contain again sub- 
subhalosJ 21 * 30 * 32 ! This is like a fractal pattern where a similar structure repeats on 
different scales which is expected to hold down to the cut-off scale of CDM. We de- 
note the different levels of subhalos as sub" halos, where n can specify a particular 
level. 

In Fig. [I] we show the dark matter mass density map of GHALO. The main 
halo and the over 10 5 numerically resolved subhalos it contains are visible as bright 
peaks. While smaller halos are falling towards a larger host halo, they get tidally 
stripped and loose mass. This lost mass is then detectable in phase space as coherent 
tidal streams that are streaming of the Lagrange points of the subhalo (see also Sec. 

a The original version of GRAFIC2 contained a bug that lead to incorrect initial conditions. Via 
Lactea II and GHALO used the corrected version of the code and are not affected by this. For 
more details consult Ref. 30. 
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Table 1. Summary table for Via Lactea I & II and GHALO. In the first 
block, basic halo properties at redshift 2 = are presented. The definitions 
of r 200b i ^200b i r v m>x i Vmax and Cy are given in the main text . The second 
block summarises some simulation properties: m p is the particle mass in the 
high resolution region, A^oob is the number of particles within r^oobi Z IC 
is the redshift of the initial conditions, e is the softening length and Lbox is 
the size of the cubic cosmological box with periodic boundary conditions. 



Simulation 


Via Lactea I 


Via Lactea II 


GHALO 


f200b 


[kpc] 


389 


402 


349 




[M ] 


1.77 x 10 12 


1.92 X 10 12 


1.27 x 10 12 


r V max 


[kpc] 


69.0 


59.8 


49.5 


Vrnax 


[km s -1 ] 


181 


201 


153 


cv 


m 


2.59 x 10 3 


4.25 x 10 3 


3.52 x 10 3 


m p 


[M ] 


2.09 x 10 4 


4.10 x 10 3 


9.94 x 10 2 


-^200b 


[1] 


8.47 x 10 7 


4.68 x 10 8 


1.27 x 10 9 


ZIC 


[1] 


48.4 


104 


57.6 


c 


[PC] 


90 


40 


61 


^Box 


[Mpc] 


90 


40 


40 



J 21 l 22 l This tidal mass loss is not complete though and the inner, very compact 
cores of these infalling halos survive as sub n halos within the host halo. 



3.2. Basic properties 

Defining the size of a CDM halo is a non-trivial task. The traditional method was 
based on the simple picture of the spherical collapse model.^This model predicts, 
that after virialisation the virial radius r v - lr encloses a region that has an overdensity 
of exactly 187r 2 (« 178) times the critical density p crit = |^ H in an Einstein-de 
Sitter cosmology. For the ACDM model the overdensity is approximately 100 with 
respect to p cr i t But CDM halos in structure formation simulations collapse in 
a different way and the region that is dynamically influenced by the halo is much 
larger than the traditional virial radius. ^ * 18 * 44 ! 

In this brief review we use the radius r2oob to characterise the size of CDM 
halos. The region within r2oob has an overdensity of 200 times the background, mean 
matter density pw The associated enclosed mass is denoted by M 2 oob = M(r2oob) 
where M(r) indicates the enclosed mass. 

Due to the shape of the mass density profiles of CDM halos (see also Sec. 4.1 1, 
they have a well defined peak in the circular velocity function 



GM{r) 

V cilc (r) = ^ — • (1) 

The peak value V^ nax of the circular velocity function V c - n - c (r) serves as proxy for 
mass and the radius ry max , where the peak velocity V max is reached, gives an inherent 
size scale of the CDM halo. Describing CDM halo properties by these two intrinsic 



b H is the Hubble parameter and G is Newton's gravity constant. 
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Fig. 1. Dark matter mass density map of GHALO. The height of the picture corresponds to 
400 kpc. The main halo, which is resolved with more than 1 billion particles, and the over 10 
numerically resolved subhalos it contains are visible as bright peaks. 

scales avoids the trouble of denning an outer boundary of halos which is of especial 
advantage when characterising sub"halos. The main properties of Via Lactea I & 
II and GHALO are summarised in Table [TJ 

4. Mass Distribution 
4.1. Mass density profile 

The simplest way of characterising CDM halos is by measuring their spherically 
averaged mass density profile. Over a large range of CDM halo masses, their mass 
density profiles can be well fit by a generalised NFW^ (GNFW) profile given by 

(r/r s p(l + r/r s y T 

This form of the profile is a special case of the af3j profile with a = 1 and (3 = 3.BSE3 
In the outer part, the GNFW profile is oc r -3 whereas in the inner part oc r -7 , r s 
is the scale radius and p s the scale density. 

The inner slope 7 was proposed to be 7 = 1, the NFW case, over a range of 
radii from 0.01 r2oob to r2oob- But it was soon realised, that CDM halos have 
some significant halo-to-halo scatter and most halos are significantly denser in their 
inner parts and have a steeper slope with 7 > 1 | 2i | 30 | 48 | 49 | 50 | 5i | 5 2j53j Thi g j s g^ju t rue 
for our highest resolution case, GHALO, which shows a steep slope of -1.4 down to 
approximately 2 x 10~ 3 r 2 oob (see Fig. [2]). 
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Fig. 2. Mass density profile (left) and slope (right) for GHALO as a function of radius. For 
GHALO r2oob = 349 kpc. The vertical line marks the convergence radius of 120 pc. 



Below this scale the mass density profile becomes significantly flatter due to the 
higher resolution and a newly proposed functional form for the mass density profile, 
the SM profile, 



(3) 



fits better The SM profile approaches the central maximum density po as r — > 0. 
We also note that the derivative dlnp/dln(l + r/R\) versus ln(l + r/R\) has a 
slope of exactly -2A. In the innermost region that is still resolved by GHALO, the 
logarithmic slope is -0.8 (see Fig. [2|.^ Such a shallow inner logarithmic slope on 
similar scales is also found in the Aquarius simulation.^ 

Locally, when measured in small volumes, the mass density can vary by orders 
of magnitudes in the outskirts of CDM halos. This is to some degree due to the 



prolate shape of the halos (see Sec. 4.3) but also because of the presence of both 
subhaloes and underdense regions or holes in the matter distribution.^ In the inner 
region, the dark matter distribution is very smooth and the contribution of subhalos 
is small™! 

Unfortunately, subhalos are still not resolved with enough particles in order to 
give a definitive answer about the shape of their mass density profile. But results 
from Via Lactea II show that they are also showing a steep inner cusp so that they 
would be similar to their host halos.^ 
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4.2. Concentration 

An intrinsic and general measure for the concentration of halos is given by the 
enclosed density within TV max in units of the critical density,'^^ 



Cy ,feJ = 2 /IW) (4) 



This definition has the advantage that it is also well defined for isolated halos as 
well as sub"halos. The traditional virial concentration 

C v ir = Tvir/Ts (5) 

depends on the detailed definition of the virial radius, the profile fitting procedure 
and has the drawback that it artificially increases due to the growth of r v i r with 
time.lSlThe usage of the inherent concentration cy circumvents these complications 
and offers a robust concentration measure of halos as long as the peak of the circular 
velocity curve is resolved. The concentrations cy for the host halos of Via Lactea I 
& II and GHALO are also given in Table [l] 

For isolated CDM halos of a given mass, objects that form earlier have a higher 
concentration. But there is a significant scatter in this relation.^ Since in a ACDM 
cosmology smaller halos form in average earlier than larger halos, they have a higher 
concentration.'2S13£l Also, isolated halos in dense environments are more concentrated 
and exhibit a larger scatter than halos in a low density environment.^ 

Generally, subhalos have higher concentrations than isolated halos.^ Subhalos 
that are located closer to the host halo center have higher concentrations than 
subhalos in the outskirts of the halo. 14121] This is due to the tidal forces that remove 
material of the infalling subhalos from the outer, low density region. Subhalos closer 
to the host halo center loose more mass due to the stronger tidal forces and more 
pericenter passages which leads to an increased concentration.l^EZ] 



4.3. Shape 

The mass distribution of isolated CDM halos has the shape of triaxial ellipsoids in 

general, though the shape is close to prolate since the lengths of the intermediate 

axis b and the minor axis c are similar E 6 I 3[| I 58 I 59 I L ess massive halos are in general 

rounder, i.e. the ratio c/a, where a is the major axis length, is larger for less massive 
halos.EHMl 

Overall, GHALO has a prolate shape as well with b/a « c/a w 0.5. The outer 
parts of GHALO, and CDM halos in general, are rounder and the axis ratios become 
larger with distance. At the halo centers the shape diverges fast and becomes more 
spherical due to numerical reasons.^ 

Subhalos have a triaxial shape in general as well but are over all more spherical 
than the host halos. Like isolated halos, subhalos tend to be less spherical in their 
central regions. Also, subhalos are tendentially slightly more spherical closer to the 
host halo center.^ 
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5. Subhalos 

5.1. Spatial distribution and orientation of subhalos 

When halos fall into a host halo and become subhalos, their mass is removed from 
outside inwards by tidal interactions. This process leads to an antibiased spatial 
distribution of subhalos with respect to the smooth background mass and subhalos 
have a more extended spatial distribution. The tidal disruption process is not able 
to completely disrupt the subhalos though and the dense cores survive. For example 
in Via Lactea I, 97% of all subhalos at redshift z — 1 survive until the present time 
and the total mass in resolved subhalos is 5 % at redshift z = 0J 14 l 17 l The volume 
occupied by subhalos increases with distance from the center of the host halo.^ 

When selecting subhalos by their current mass, they have a number density 
distribution at redshift z — of 

™ S ub,M(V) oc p{r) r (6) 

where p is the smooth DM density. ^ 4 * 32 * 60 * The subhalo number density is indepen- 
dent of the selection mass threshold and the subhalo mass, i.e. subhalos in different 
mass bins have the same spatial distribution when properly normalised.* 32 * 60 * ^ Vmax 
is less affected by tidal mass removal and therefore V max selected subhalos follow 
closer the smooth DM distribution at redshift z — with 

"■sub,v max (V) oc p(r) M(r) (7) 

where M(r) denotes the enclosed mass.E^- 21 * 32 * 

During their orbit within the host halos, the major axis of subhalos tends to be 
aligned with the direction toward the host halo center. This is effect is most likely 
due to tidal interactions as well in the sense that tides stretch the subhalo shape 
along the direction towards the center. For subhalos further out (i.e. distances larger 
than 3 ?~2oob), this alignment is not observed any more.^ 

5.2. Sub n halo abundance 

The cumulative subhalo abundance within r 2 oob when selected by velocity (velocity 
function) at redshift z = can be described by 

N{> U oc (v^S«_ (8) 

\ * max, host / 

with ay = 3 and where V, na x,host is the peak velocity of the host halo.*^ 2 - 1 * 32 * 49 * For 
Via Lactea II the proportionality factor is 0.036 and \4nax,hoBt = 201 km _1 (see 
also Table [TJ .I3H The velocity function is independent of the selection radius, i.e. 
in the case of Via Lactea II the velocity functions of subhalos selected within 100 
kpc or 50 kpc have the same shape as the velocity function for subhalos selected 
within r2oob, only the normalisation changes.* 2 ^ The mean sub 2 halo abundance is 
consistent with a scaled down version of the main halo within the current resolution 
of simulations.^ 
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When selecting by mass, the cumulative subhalo abundance (mass function) at 
redshift z — is well described by 

f M \ _QJU 
N(> M) <x _ (9) 

\M 200b, host / 

with a M = 0.9 . . . l.O.' 13 ' 14 ' 32 ' 48 ' 49 ' 

6. Velocity Space Structure 

6.1. Anisotropy parameter 

A parameter to describe the velocity space structure of CDM halos is the anisotropy 
parameter 

/3=i- ±4* (io) 

1 CT rad 

where ot an is the tangential velocity dispersion and cr ra d the radial velocity dis- 
persion. A value of (3 « means isotropic, (3 > respectively (3 < denotes a 
radial respectively tangential anisotropy. For relaxed CDM halos (3(r), measured in 
spherical shells, changes from nearly isotropic in the center to slightly radial in the 
outskirts of the halo.' 22 ' 50 ' 60 ' 62 ' 63 ' 64 ' 65 ' 

When (3 is measured in small, local volumes, then the local value can greatly 
deviate from the spherically averaged value. In the inner halo, regions along the 
long shape axis particles are preferentially on radial orbits whereas along the short 
shape axis particles are on tangential orbits. This effect disappears at larger radii.^1 

6.2. Velocity dispersion ellipsoid 

As a consequence of the tensor virial theorem, the global velocity dispersion ellipsoid 
is aligned with the shape ellipsoid.^ This is also true for small local volumes in the 
inner parts of halos, where the local velocity ellipsoid is aligned with the global 
shape as well. In the outskirts this local alignment is not observed any more.^ 

The shape of the global velocity dispersion ellipsoid (and as a consequence the 
potential) is significantly rounder than the mass distribution.^ When measured 
again in small local volumes, the shape depends on the location within the halo.^ 

6.3. Velocity distribution function 

From theoretical models that describe CDM halos with smooth phase space distri- 
bution functions, it is expected that the local velocity distribution has the shape 
of a normal distribution around regions where the mass density profile has a slope 
of -2, i.e p oc r~ 2 . In the inner regions, where the mass profile is shallower, the 
velocity distribution is more peaked and in the outer region, where the mass profile 
is steeper, the velocity distribution is broader J 46 * 66 * 67 ! 
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r = 8 kpc r = 400 kpc 




Fig. 3. Velocity structure of Via Lactea II. Positions (top row) of particles within small volumes 
at 8 kpc (left) and 400 kpc (right) from the Via Lactea II center as well as the corresponding 
location of the particles in the local velocity space (bottom row). The colour in both cases is given 
by the magnitude of the velocity vector. 

But a smooth velocity distribution is only a good description of the velocity 
space of CDM halos in the central region. In the outer parts, the local matter stems 
from several tidal streams of subhalos that lost material while passing. This is best 
illustrated in Fig. [3] This figure shows the positions and velocity vectors (top row) 
for every particle within local volumes centered at 8 kpc (left) and 400 kpc (right) 
from the center of Via Lactea II. We also plot the location of the particles in the local 
velocity space (bottom row) . The colour (and length in position space) encodes the 
magnitude of the velocity vector. The big white arrow in the center of the position 
space plots points towards the galactic center whereas the black cube in the center 
of the local velocity space plots marks the origin. For the velocity space plots, we 
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Fig. 4. Central sphere of radius 50 kpc where we calculated the position space density (left) and 
the true phase space density (right) with EnBiD. Approximately 2000 peaks from the subhaloes 
within that radius are visible in both pictures but the contrast of subhaloes is higher in the phase 
space density picture. In the phase space picture are also large scale dark matter streams visible 
which were formed by tidal mass loss of infalling subhaloes. These streams are not visible in the 
traditional position space density picture. 

split the velocity vector field in a radial- (x-axis), ip- (y-axis) and a -^-component 
(z-axis). 

These small volumes have been selected to be free of subhalos. Therefore, one 
would not expect to see a clumpy structure in the velocity space of these volumes. 
At 8 kpc no velocity space structure is apparent. The orientation of the velocity 
vectors appears random and the actual velocity distribution can be well fit by 
a multivariate normal distribution. But at 400 kpc, there is evidence for locally 
coherent motion, visible as groups of vectors with the same colour pointing in the 
same direction. This is even clearer in the velocity space plot: there are only a 
bit more than a dozen clumps in velocity space and no smooth component at all. 
Obviously a smooth multivariate normal distribution would not provide a good fit. 
This is evidence that the outer part of the halo is built up from a collection of 
large scale streams from the tidal disruption of infalling subhaloes. But so too is the 
smooth component in center, in that it also likely consists of an overlap of many, 
many streams. 68 69 The central limit theorem then guarantees that the resulting 
distribution in the center closely resembles something like a multivariate normal 
distribution^ or more general distributions .1221121 

7. Phase Space Structure 

In Fig. [3] we show the central region of Via Lactea II: the sphere has a radius of 50 
kpc. We plot the position space density (left) and the true six dimensional phase 
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Fig. 5. Phase space density profile for GHALO. The Q(r) /(2n) 3 ' 2 profile follows a power law 
with slope -1.84 (dotted line). The mean of the true phase space density calculated with EnBiD is 
dominated by the contribution from subhalos beyond 1.8 kpc. The approximate power law relation 
can be recovered by removing high phase space density contributions from subhalos and streams. 
The true phase space density relation is steeper and has a slope of ca. -2. The median of the true 
phase space density is not affected by the contribution from substructure. The resolution limit of 
GHALO is 120 pc and is marked with a vertical line. This figure was adopted from Ref. 30. 



space density (right) calculated with EnBiLj^J^D We only show the top five to six 
orders of magnitude in position and phase space density, in order not to overload the 
two pictures. Many cold streams are clearly visible in the central region. Although 
these cold streams only contribute a few percent to the local mass density, their 
velocity dispersion is just a few km s _1 , resulting in a very high phase space den- 
sity for these particles. These streams are not visible in traditional density picture 
(see Fig. [IJ and can only be revealed by visualising the true phase space density. 
Via Lactea II seems to be the first structure formation simulation with sufficient 
resolution to reveal these streams in phase space, since previous lower resolution 
runs did not show such phase space features. Also subhaloes are better visible due 
to their higher contrast in phase space density. Approximately 2000 peaks are seen 
in the central 50 kpc sphere phase space density image. The high contrast makes 
the phase space density the ideal method of identifying subhaloes. 

An often used method of describing the structure of phase space is by the pseudo 
phase space density Q = p/u 3 where p is the matter density and a is the velocity 
dispersion. When Q is measured in spherical shells, then it follows a power law 

Q(r)cxr" a « (11) 

with an = 1.8... i.95.™mE3lHI75] 



c http : //sourcef orge .net/pro jects/enbid/ 
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In Fig. [5] we plot the quantity Q(r)/(27r) 3 / 2 which follows a nice power law 
with a slope of ctQ = 1.84 for GHALO^It is interesting to compare this to the 
spherically averaged true phase space density again calculated with EnBiD. When 
taking the mean true phase space density in spherical shells, we see that subhalos 
start to dominate the profile from ca. 1.8 kpc on where the innermost subhalo 
in GHALO is located. In the inner part, the mean EnBiD value roughly shows 
a power law behaviour but with a steeper slope of -2. Subhalos have in general 
a higher central phase space density than the host halos.^^1 By cutting out the 
contribution of particles with / > 100 M Q kpc~ 3 km~ 3 s 3 we can remove the 
contribution from subhalos. When removing also streams by neglecting particles 
with / > 0.4 M© kpc~ 3 km~ 3 s 3 the approximate power law behaviour can be 
extended out to ca. 40 kpc. Also, the median of the true phase space density is not 
affected by the contribution form substructure (see Fig.[5|. 

The power law behaviour of the pseudo phase space density Q is still enigmatic 
and likely reflects the distribution of entropy K = Q~ 2 ^ 3 = a 2 /p 2 ^ 3 , which dark 
matter acquires as it is accreted.^ Therefore, one should use Q with caution and 
give preference to the true phase space density when studying the phase space 
structure. 

8. Conclusions and Final Remarks 

In this brief review we discuss recent results from the high resolution structure for- 
mation simulations Via Lactea I & II and GHALO which study the formation of 
Milky Way sized objects. The resulting CDM halos show a wealth of structure: sev- 
eral levels of sub ra halos and streams in phase space in regions that where smooth in 
previous simulations. We present quantitative descriptions of the mass and subhalo 
distribution as well as a detailed insight into the velocity and phase space structure. 

Unfortunately, the content of this brief review is limited and we had to skip 
results on implications for indirect and direct dark matter detection experiments. 
Please consult Refs. 13, 23, 25, 26, 28, 29 for more results in that area. For a 
complementary review on the structure and evolution of CDM halos consult Ref. 
78. 

Images, movies and data from Via Lactea I & II and GHALO are available at 
http: //www.ucolick . org/^diemand/vl/ and http : //www. itp .uzh. ch/ghalo/. 
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